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Abstract 

We study uncharged Rindler hydrodynamics at second order in the derivative expansion. The 



CN . equation of state of the theory is given by a vanishing equilibrium energy density. We derive 

relations among the transport coefficients by employing two frameworks. First, by the requirement 

j-v of having an entropy current with a non-negative divergence, second by studying the thermal 

(D ■ 

r~| ' partition function on stationary backgrounds. The relations derived by these two methods are 



consistent with each other. However, we find that the entropy current yields stronger constraints 
than the thermal partition function. We verify the results by studying explicit examples in fiat 



> 

in 

^y-N ■ and curved space-time geometries 
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I. INTRODUCTION AND SUMMARY 

There has been much interest in recent years in the holographic relation between the 
hydrodynamics of quantum field theories defined on [d + l)-dimensional space-time and 
deformations of black hole geometry in one higher space dimension. Aspects of this relation 
have been studied for a large class of field theories in diverse dimensions including conformal 
and non-confomal theories, relativistic as well as non-relativistic ones, anomalous theories 
and theories with spontaneously broken symmetries. 

One intriguing example, which will be the subject of this paper, is the Rindler hydrody- 
namics, that is the hydrodynamics induced on a codimension one hypersurface in Rindler 







geometry [m]^- The Rindler geometry is a solution to the vacuum Einstein equations, and 
correspondingly the equilibrium energy density in the hypersurface hydrodynamics vanishes. 
Despite this unconventional equation of state, the hydrodynamic description is well defined. 
Viewing hydrodynamics as an effective description of quantum field theory at large length 
scales and at local thermal equilibrium, one may hope that understanding Rindler hydro- 
dynamics can shed light on how Holography works in asymptotically fiat space-times. 

In the hydrodynamic description one upgrades the local charge densities to local cur- 
rents. In this framework, the entropy current in hydrodynamics is an upgrade of the local 
entropy density to a local current [5,]. The entropy current has not been given yet a mi- 
croscopic origin and it is not a unique object. It is constructed phenomenologically order 
by order in the hydrodynamic derivative expansion. One can generalize the second law 
of thermodynamics, by requiring that the entropy current has a non-negative divergence 
order by order in this expansion. The requirement imposes constraints on the transport 
coefficients of the hydrodynamic description, which in turn impose constraints on the un- 
derlying microscopic quantum field theory. This has been applied, for instance, in the study 
of anomalous hydrodynamics gLZJ, uncharged hydrodynamics |8!], and in the analysis of 
holographic hydrodynamics p. Il0l|. 

There are two types of such constraints. One type gives equality constraints that relate 
different transport coefficients. The second type gives inequalities, such as the requirement 
that the shear and bulk viscosities at the first viscous order are non-negative. In order to 
fully exploit these constraints one studies them in the presence background fields. In the 
case of uncharged hydrodynamics the background field is a curved metric. 



A second method to derive the first type of constraints employs the thermal partition 



function on a curved background with a timelike Killing vector ll|, |l2|. In this framework 
the constraints are related to the underlying gauge and diffeomorphism symmetries. 

The aim of this paper is to study the first type of these constraints in the case of Rindler 
hydrodynamics. We will employ the two methods discussed above and compare the results. 
The paper is organized as follows. In section 2 we will classify the independent fluid data, 
i.e. modulo the conservation laws of the stress-energy tensor in ideal hydrodynamics. We 
will use as variables the fluid velocity and pressure as well as the curved background data. In 
section 3 we will construct the entropy current up to second order in the derivative expansion 
and its divergence. By requiring the latter to be non-negative we will derive constraints on 
the transport coefficients. We will check the constraints in two cases: flat hypersurface 
(with and without a bulk Gauss-Bonnet term) and curved hypersurface in Einstein gravity. 
In section 4 we will constrict the thermal partition function and derive the constraints on 
the transport coefficients. We will compare the two and see that the latter are somewhat 
weaker. Some details of the calculation are outlined in the appendices. 

II. THE FLUID DATA 

We would like study the (d+l)-dimensional Rindler hydrodynamics in a general curved 
background up to second order in the derivative expansion. The fluid variables that we 
will use are the fluid velocity m^ normalized as m^m'^ = — 1, and the fluid pressure p. The 
curvature data will be given by the Riemann tensor R^p„. In Rindler hydrodynamics we have 
a vanishing equilibrium energy density eo = 0. We will start by classifying the independent 
fluid and curvature data, that is, modulo the ideal fluid hydrodynamic equations and the 
curvature identities. The classiflcation of fluid and curvature data when eo 7^ has been 



carried out 



in [81. 



One cannot simply use this classiflcation by setting eo = 0, since the ideal 
fluid hydrodynamic equations are different. We will use, however, the same classiflcation 
scheme according to the transformation properties under the local SO{d) symmetry group 



a 



that leaves the (d+l)-velocity u^(x) invariant, as employed in 

The Rindler hydrodynamics stress-energy tensor at the ideal level is 



ri? = PP,U , (1) 



and the conservation equations VJT'"^ = projected on Ua and on P^ are: 



e = v^u^ = , 



U^'VaU" + P^"VMp = . 



(2) 
(3) 



pt^'^ = gt^'^ -\- u^u^ is the projector on the vector space perpendicular to u^. The fact that 
the fluid expansion vanishes will be of importance. 

In the following tables we classify all the possible terms after imposing the ideal fluid 
conservation equations. The independent fluid data at first order is given in Table [H 





Independent data 


Scalars (1) 


D\np 


Vectors (1) 


a^' = {u.V)u'' 


Tensors (2) 





TABLE I: Fluid data at first order in derivatives 

D = u'^V^, a^ is the fluid acceleration, a'^" and u'^'^ are the fluid shear and rotation 
tensors, respectively. We define A(^^By) = }^{A^By + A^B^) and A^^Bi^-] = ^[A^B^ — A^B^] . 
At this order there is no curvature data. The independent fluid data at second order is 
derivatives in given in Table HIl The composite expressions of fluid data that we will need 





Independent data 


Scalars (1) 


D{Dlnp) 


Vectors (2) 




Tensors (1) 


P/^^P'^^VaV^lnp 



TABLE II: Fluid data at second order in derivatives 



for the analysis are given in Table lllli We define 

A^ ^ _ papfS fAafS + Apa 



9aP 



PP^A 



pa 



d 



(4) 



Scalars (4) 


{D\np)\ o2, a;2, a^ 


Vectors (3) 


a^'D\np, ayoo^"" , aua^"' 


Tensors (5) 


Dlnpa^„, 


^f^f^au), '^^^^au)i ^^^^au), 


0-{nO-v) 



TABLE III: Composite expressions in fluid data at second order 

The classification of the independent curvature data is the same as in |8| and for completeness 
we give it in IIVI 





Independent data 


Scalars (2) 


Roo = u^u^Rfjiu 
= u'^u'^R'^^^, 


Vectors(l) 


P^'^R^puP 


Tensors(2) 


where K^y = u^ul^R^avp 



TABLE IV: Independent type curvature data at second order 



Consider next the the most general form of the stress-energy tensor of Rindler hydrody- 
namics in curved space. We fix the ambiguities at the derivative order by requiring that 







y 



P^Tl-\^ 







(5) 



and the pressure receives no derivative corrections. With this choice, the stress-energy tensor 
up to second order in derivatives takes the form 



+ Al Cr^'^aaiy + A2 Cr(^°WQiy) + A3 Uf^°'Uau + A4 a^fli- 



(6) 



+ UuU 



fl"'U 



dia"^ - t/aw^ + ds{D\npy + diD^Dlnp) + rfga^ + eiR + Cgi?! 



■00 



where 



K"" = R'''"'^uc,u^, R"" = R'^'^'^gap • (7) 



There are two transport coefficients at ffist order: rj is the shear viscosity and C is a con- 
tribution to the energy density. There are fifteen transport coefficients at second order. 
This is the same number as for a generic non-conformal hydrodynamics with non- vanishing 
equihbrium energy density. However, some of the transport coefficients muhiply different 
expressions in the fluid variables. In particular, seven of these terms can be viewed as correc- 
tions to the energy density. As we saw in the previous section, the fluid expansion O = V^m'^ 
vanishes and the scalar Dlnp appears instead in the stress-energy tensor. Note, however, 
that unlike the case where the equilibrium energy density is nonzero, here this scalar is not 
related to 9. 

III. THE ENTROPY CURRENT 

In this section we will construct the most general entropy current that exhibits vanishing 
equilibrium energy density. We will construct it up to second order in derivatives and derive 
the constraints on its form and on the transport coefficients in the stress-energy tensor. 
These follow from the requirement that its divergence is non-negative. 

A. The entropy current at second order 



B 



In [81], the most general entropy current when the equilibrium energy density is nonzero 
has been constructed. In the following we will use a similar procedure and outline the 
differences. As above, we will use as fluid variables the normalized fluid velocity m^ and the 
pressure p, and replace the fluid expansion 9 by the scalar Dlnp. 



The entropy current J^ has a derivative expansion 

where Jf^-, = su^ and J^n = 0. At second order in derivatives we get that the most general 
entropy current depends on thirteen parameters, which correspond to the six independent 
vectors (three of them are composite), and seven independent scalars (multiphed by u^) at 
second order, as outhned in the tables. It is convenient to parametrize the entropy current 
at the second order as follows: 

Jf2) =V, [Ai(M^V"lnp - M^V'lnp)] + V, (Aslnp u^'"') 

+ A3 (r^" - i(7^"i? ) u, + [A^D{D\np) + A^R + Agi^oo] u^" 

\ ^ J (9) 

+ {Biu'^ + B2{Dhipf + Esa^)^^ + 54 [a^M'^ - 2(Dlnp)V^lnp] 
+ [DlnpV^^s - P"^(V^u'^)(Va55)] + EgZ^lnpa^ + fiya^a'^'^ . 

The divergence of the entropy current up to second order reads ( TT2|) : 

V^J^ = V^(sn'^) + V^Jf2) 

= --a,,W'' - -De + V^ Jf2) . (10) 

where the energy density at second order e = T^u u'^u'^ is 

e = rfio-^ - rfsw^ + c?3(L'lnp)^ + (i4D^lnp + d^a^ + eiR + Cai^oo ■ (H) 

Note, that the entropy density s is constant up to second order since it is related to the 
energy density by the thermodynamic relation de = Tds. 

B. Constraints on the entropy current and the transport coefficients 

In order to constrain the entropy current and to find relations among the transport 
coefficients in the stress-energy tensor, we have to analyze the terms that can appear in the 
divergence of the entropy current. We omit in the entropy divergence terms that cannot 
be completed to a positive expression, that is to a complete square. We will encounter two 
cases. In the first case we find terms that cannot appear in the second order divergence 
of the entropy current, for instance, [Dlnp)"^. We will therefore set to zero the coefficients 

8 



of all the third order terms in the entropy current divergence that take the form of second 
order scalar times Dlnp. In the second case we identify terms that cannot appear at the 
fourth order divergence such as {R^vY and set to zero the coefficients of third order terms 
that require this term for a complete square, e.g. a^'^R^y. For this we will have to construct 
the most general entropy current to third order in derivatives and take its divergence. Note, 
that in order to simplify the analysis we will use in the calculation P^'^P^^Daap instead of 
pA^^p^/^V^V/^lnp. 

Carrying out this analysis, we find the most general entropy current with positive diver- 
gence : 



+ A3 (r^'' - \g^''R\ uu + A,D{D\np) - YJ^R + A,R 



-e-rtoo 



W 



-^)u^' 



( dA, v'd'A.l vdA,\ (12) 

+ (a^ + p'^^^ [DhipvnYip - p'^^iy pu^'){y Mp)] 

( , dA-x ^d^A-i , dAa\ , 

\ dp dp'^ dp J 

For a detailed computation of its divergence see Appendix [XI 

The relations among the stress-energy tensor transport coefficients that we get are: 

C = (13) 

dKi 

K2=T-Ki+p—- (14) 

dp 

pdKi 
pe, = -1 - 2^ (15) 

A4 = p{2ei - 62) + P-i-i'^ + Ki- K2) -2t - Ki + K2 (16) 

7) d 
Q = pd2+ pei + pe2 - 2ki - k,2 + h + t^-(2/«i + 1^2- A3) (17) 

4 dp 

pd5=pei-pe2 + p^ -^{261-62) (18) 

dp 

We see that the energy density can only be corrected at second order. 



C. Flat space-time background 

In the particular case of hydrodynamics in flat space-time, we set to zero the curvature 
terms in the entropy current and in the stress-energy tensor. The relations among the 
transport coefficients of the stress-energy tensor that follow from the requirement of non- 
negative divergence of the entropy current are 

C = (20) 

6 dr d'^T 6dci , , 

p dp dp^ p 

(22) 



4d2 

p 

20?3 


^dd2 ^ 8A3 
dp p^ 

d fdA 


SdAs 
p dp 


8A4 ^ lrfA4 

p^ p dp 


16r 
p2 


p' 


dp\p) 









The coefficients of the entropy current are restricted by 

d_ 
dp 



57 = . (25) 

The relations fl20tl25p can be checked against two examples of Rindler hydrodynamics derived 
via the gravitational solutions in one higher space dimension. The first example is the fiuid 
hydrodynamics obtained in the framework of Einstein's gravity, i.e. from a gravitational 
background obeying the Einstein equations for an empty universe Rab = 0. This example 
was studied in |2|, yj, and its results were: 

--Be = 2sp-^ B; = 3sp-\ p^ = 2sp-\ Bj = 

dp 

di = d^ = d = €2 = 0, 

= -4, Ai = -2 

(26) 

s = 47r is the entropy density in equilibrium in this case. 

The second example is an empty universe in Gauss-Bonnet gravity. This case has been 

worked out in |^, its results are the same as in the previous example with the following 

10 



Ai 


= '^sp-\ 


A2-- 


= -^3 


Bi 


1 -2 
= -^sp , 


B2 


= B,= 


di-- 


= -2p-\ 


d2 = 


--d,= 


T = 


Aq = A2 


= A3 


= -A4 


v = 


1, C = 


. 





modifications: 



A2 = -A3 = sp 2(1 + 4a/), A5 = --sp 2(1 + 4V), Aq = sp 2(1 - 4ap2), 
5i = ^sp-2(l + 4V), fi2 = 56 = 2sp-2(i_2a/), 53 = sp-2(3 + 2ap2), 

p— ^ = 2sp~2n _ 2^^ 
dp 

A3 = -4(1 + 3ap), Ai = -2(1 + 2ap) , (27) 



where a is the coefficient of the Gauss-Bonnet term. In both cases the relations (I20ti25p are 
satisfied. 

D. Curved space-time background 

In this subsection we will construct a solution in Einstein gravity that describes holograph- 
ically Rindler hydrodjTiamics in a curved background. The additional space coordinates will 
be denoted by r and the hydrodynamics is defined on the co-dimension one hypersurface 







r = Tc, We will follow the procedure in [3| where a fiat metric ri^.^ was chosen on r = Tc, and 
we will replace it by a curved metric g^^. The details of the construction are presented in 
appendix [Bl We outline them in the following. We need to construct the bulk metric up to 
second order in derivative expansion. We do that by solving the bulk equations Rab = 0, 



requiring regularity of the so 
The difference compared to 



ution at r = and fixing the metric on the hypersurface r = Tc- 
3| appears at second order in the derivative expansion, because 

curvature effects are of second order or higher. The second order correction to the bulk 

metric reads 

^curvGf] = 2u^u,ir - r,)R - 2p\r - r,)\^Pi;^R^^ux + 2(r - r,){R^, + K^,) , (28) 

where i?, R^^ and K^y refer to the curvature of the hyper surface metric. The entropy (area) 
current J^ = -^^y/Gl^, where Gn is the Newton constant, is affected by the corrections of 
the metric determinant and the normal to the event horizon: 

6,urvG^^'^ = --{R + Rm) , (29) 

P 

5curvt^2) = -^P^'RS ■ (30) 



11 



The final result is the same as (126|) with the following modifications: 

As = As = -2sp"^ ei = -2p-\ Ki = K2 = -2 . (31) 

It is easy to see that this result satisfies the constraint equations: fll3H19l) . 



IV. THERMAL PARTITION FUNCTION 



In this section we will consider the Rindler thermal partition function on curved back- 



grounds with a time-like killing vector Uj, |l2|. We will use the framework employed in 



ll| to find relations among the transport coefficients for uncharged non-conformal hydro- 
dynamics. As above, since the equilibrium energy density vanishes in our case, we cannot 



simply borrow the results of ll|] that assume its non-vanishing. Indeed, while the relations 



among the transport coefficients derived in ll| are the same as those obtained from the 



non- negativity of the entropy current divergence [8| , this will not be the case here. 

The equilibrium backgrounds posses a time-like killing vector and can be written in the 
Kaluze-Klein form 

dsj^j^ = G^^dxf'dx" = -e^'^^^) {dt + ai{x)dx')^ + gij{x)dx'dx^ , (32) 

where i = l...d. 

The fluid variables take the form 



4 = e-'^(-^)(l,0,. 


•,o), 


P = Poe-''^'^ 




eo = . 





(33) 
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A. The stress-energy tensor 

Evaluating the derivative part of the stress-energy tensor ([6]) on the curved equihbrium 
background we get 

11% = a, (b!J - ^A + a2 (v W - V'agA + «3 (vVVV - ^^^gA 

+ a, (f'fj + ^gA e'- + g'^ Lr, + hV'a + hiVaf + w\feA (34) 
where, 



Jij (Ji<Xj OjCLi, J JijJ 

blP = ^f^l, hp = K2- Ki-T 

hp = 2*^'^2 - 1^1 + A4) 

hp = t(2/«i + K2- A3) , aip = Ki 

a2P = K2- Ki-T, a-ip = /t2 - Ki + A4, a^p = -t(2ki + ^2) + -A3 . (35) 

K^JjRg are Ricci tensor and Ricci scalar respectively, calculated from the spatial metric 
gij{x) in ( l32l) . The corrections of the fluid variables to second order are: 

u'' = bou^ + Vie'^V^aP' + V2e''V,r\ 
60 = 1 - {v^e^'^a.V.ar + V2e^'' a,V , P') 
P = Poe"" + hRg + t2VV + t3(Va)2 + ^47/^62^^ 



4' 
e = nRg + r2 VV + rs{Vaf + n^fe^^ , 

where 

n = ei, r2 = 62-61, r3 = 62 - 61 + (is, r4 = -(26i + 62 + 6^2) . (36) 

The second source of corrections arises from inserting the velocity and pressue corrections 
into the zeroth order stress-energy tensor. We find that the modification of the stress-energy 
tensor due to these corrections is given by 

"2„2cr 



T' = g'^ ihRg + t2VV + t3(Va)^ + U-pe 
Too = e'^ (ni?, + r2 VV + r,{Var + n\fe'^^ (37) 

T^ = -Po {vie'-V.aP' + t;26'^V,/^'*) 

13 



Adding ([M]) to ([37D we get 

T'' = 9'' (hRg + t2 VV + h{Waf + U-J^eA + Ii% 

Too = e'" (nR, + r^V'a + r^iVaf + n\fe^'''\ (38) 

where 11'^^ was listed in (l34l) . 

B. The partition function 

The thermal partition function takes the form 

W = \ogZ = -]^j d'^x^, [S,Rg + pIS^U.P + S^iVaf] , (39) 

where 5*1, 5*2, S^, are three arbitrary function of a. 

In order to evaluate the stress-energy tensor from the partition function, we vary the 
partition function 

where we integrated over the time coordinate t which is compactified on a circle with radius 
— . One obtains the stress-energy tensor by varying with respect to the metric, 

^ Poe^" 6W ^■_ Po SW_ 

T- = ^^^,V-|^. (41) 

The stress-energy tensor derived from the thermal partition functions (139|) reads 

sT^ = pS, (i?y - ^Rgg'^) + 2pIpS, {f'f,, - \fg'^) + p{S, - S'i) (VVVV 
- ^(Va)V^) -pS[{V'Wa - g^^V^a) +]^pS';{Vafg^^ 

sToo = f^{S[Rg+plSU^-S',{Vaf-2S,V'a)) 
sT^ = 2plp{S',V,ar + S2V,r) , (42) 

where ' denotes derivative with respect to a. 

14 



Comparing the equations (l38l) and (H2l) . one can express the transport coefficients in 
terms of the three coefficients S appearing in fl39l) . We find 

sai = pSi, 502 = -pS'i sa^ = 2plpS2, sa^ = p{S3 - S") 

1 , (43) 

sbi = -pti, sb2 = -pt2, sh = -pt3 + ^pS'l, shi = -pti. 

One can ehminate the coefficients S's from above set of relations which gives a relation among 
transport coefficients. These relations are the same as those we got from the positivity of 
the entropy current in the previous section. However, we find now only five (I14til8p of the 
seven relations, the relations ( 1T5|) and ( IT^ are not revealed by this method. 
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Appendix A: The entropy current divergence 



Here we present the full computation of the entropy current divergence (ITO|) : 

p p 

p— p— h /-D4 - p— h -De — 3— -t^lnpitoo 

dp dp dp dp J 



+ 



\ dp 2 dp dp J 

+ fp^ + 25, _ 2p^ - ^ - ^"l D^lnpDlnp 
\ dp dp p dp J 

f dBi ^d^B^ dB^ dBe 2d5 dd^ \ 2 ni 

\ dp dp^ dp dp p dp 

+ (A,-^]DHnp+(-p^^' 



P J \ dp 



Be + ^) a^V,(Dlnp) 



(„ „ nd B^ dBj 2d] 2d'!, „, \ o,, , , 

+2^3 - 2^4 - p'-— ^ - p—^ ^ + — - p-^K a^^y ^InpW p\np 
dp^ dp p p J t^ t- 

+ [A,~''AdR+ (a^ - ^) DRoo 

p J V pJ 

, dBi dBr. ddo 2(^2 \ 9 , 

+ p—^ - 2Bi + 2^4 - 2p—^ + Be + -^ -) co^Dlnp 

dp dp dp p J 

, dB-i ^ ^ dB^ ^ ^ ddi 2di 9 , \ 9 ^, 

+ p—^ - 2B3 + 2^4 - 2p—^ + Be-Br--^ + ^- p"^Xo cr^Dlnp 
dp dp dp p J 

-p-'\, - 2Bs + ^) a^^a^.al 
P J 

+ ( -p-\, - 253 + ^) a'^'K^x 

4^2 , 0, 2di\ y^ 



2B3 + 4B, + — + p-'X3 a'^^u^u.u 

P P ) 

^d B^ dB^ dB^ dds\ 3 

P ~rY' + P~j P~r~ + T~ \D\np) 

dp"^ dp dp dp J 

- (2B; + p^ + Bj-'^+ p-'f\ a'^'^V.VMp - S70, V^a^^ (Al) 

Appendix B: Curved background solution 

In this appendix, we explain briefly how to get the results of the stress tensor and entropy 
current for a curved hypersurface metric. We begin by solving the Einstein equations: 

16 



R-AB = 0, then boosting the solution with the velocity vector -u^. Now we promote the 
parameters of the solution to be depended on space-time coordinates of the cut off hyper 
surface located on r = r^.- Obtaining the solution: 

ds'^ = -(1 +p^(a;^)(r - rc))u^{x^)uu{x^)dx''dx'' - 2p{x^)u^{x^)dx^dr + P^^{x^)dx''dx'' . 

(Bl) 

where we choose the gauge choice to be: 

grr = 0; gr^^ = -pu^ . (B2) 

Since we promote the parameters to be dependent on the coordinates of space-time, the 
metric does not satisfy the Einstein equations. We will solve Einstein equations order by 
order in the derivative expansion. The Einstein equations at the n*'^ order are: 

^Rab + R?B = , (B3) 

where R^b i^ ^^e n*'^ order Ricci tensor calculated from the (n-1) order metric, and (^R^b 
is calculated from the correction of the metric to the n**^ order. The final result is 



5Rlr^ = -Id^iP'^'Sg^xJ) (B4) 

5RS = \p^AiP"''S9^xJ) + y'd'Au%^;^) (B5) 



„ \"fj,-i\- -JVAI • "iy-i\- - J ^A J J j r)-i\'J^V/ n 



flf 



- lu,uMu\-Sg^;'J) + ^^u.uMP'^'SgxJ)- (B6) 

where $ = 1 + p^(r — re). 

This solution needs to be consistent with the following boundary conditions: (i) no singu- 
larity at r = and (ii) a given curved induced metric on r = Vc- The second implies that 
all the n > 1 order corrections vanish at r = r^. Projecting into components normal and 
transverse to m^ we find 

n^P;^9t} = (V2)(l - r/r,)V;(")(x) - 2p j" dr' /" rfr'T^^R^^ , (B8) 

J r J r' 

u\''5g^;^} = (1 - r/r,)A(")(x) + p T dr' T dr" (pP^"Ri"j - p-^^n^^^) - 2R; 

Jr Jr' ^ 



(B9) 
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where V^ {V^ u^ = 0) and A*^"^ are free, undetermined functions at this stage. To deter- 
mine these functions we calculate the Brown- York stress-energy tensor on the hyper surface 
located at r = r^. We impose the frame condition that the pressure does not receive any 
derivative corrections. This condition eliminates the derivative parts of the term propor- 
tional to g^t, in the Brown- York stress-energy tensor and fixes A^"^' to a non-zero value. Note 
however that the stress tensor now has a term proportional to u^u^. Hence at this point 

(2) 

we see the traditional Landau frame T^^u^ = will be inconsistent. To fix V^ we instead 
require T^^^^P^^u^ = 0. 

With these values fixed, the stress-energy tensor can be put in a more conventional form: 

- 4p-ip;P,^V„V^lnp - Ap-'a^,D\np + Ap-' Df^\npD^\np - 2p-^P^P^R^fi - 2p-^K^, 

(BIO) 

Due to the u^^Uy term, one can see that the energy density is no longer zero. It is corrected 
at the second order: 

e = T^yu^v" = --a^ua^" - -R . (Bll) 

p p 

Since metric solution is no longer stationary, the event horizon is dynamical and its location 
varies in time and space. In order to find rh{x^) we will need to solve the following equation 
in the derivative expansion 

g^^dA^r - rh{x^))dB{r - rh{x^)) = . (B12) 

Using our previous results for the metric, it is straightforward to show at second order 

12 S 1 8 

rh = rc ^ + -^u^'V^lnp - — ^0-^/30^"^ - — -w^/^w"'^ -D\npD\np 

pz p.i 2p'^ 2p^ p^ 

+ \D^nnpD;:lnp + ^D{Dlnp) -^{R- Roo) • (B13) 

Now, it is easy to construct the normal to the event horizon by i^ = g'^^dBir — Vh) and the 
metric determinant. We find the entropy current 

J^ = ^ ("l - 4 (^«/3^"^ - l^a^^"^ + 2P"''V«V/3lnp - 2P"^V JnpV^lnp + R + 2Roo^ ^ 

(B14) 
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